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I. INTRODUCTION 

The phase structure of Quantum Chromodynamics (QCD) 
is the subject of very active experimental and theoretical 
research. A crucial question is how to detect the formation 
of the quark-gluon plasma in heavy-ion collisions at ultrarel- 
ativistic energies. Vector mesons play a very important role 
in this context because they provide promising evidence for 
both deconfinement and chiral symmetry restoration. While 
the suppression of heavy quarkonium could be a signature 
of deconfinement [1], in-medium modifications of light vec¬ 
tor mesons may signal chiral symmetry restoration [2] . The 
latter manifest themselves in low-mass dilepton data from 
heavy-ion collisions [3-5]. Dileptons escape the fireball 
essentially without interaction and couple directly to light 
vector mesons such as the p. Thus, dilepton spectra show 
prominent vector meson peaks which allow for the investi¬ 
gation of in-medium modifications of these mesons [6] . A 
connection between the modifications of vector mesons and 
chiral symmetry restoration in a hot and/or dense medium 
can be established e.g. by considering the scaling of the p 
mass with temperature [7-9] or the melting of the p reso¬ 
nance [10, 11]. This connection is based on the fact that 
chiral symmetry restoration implies the degeneration of chi¬ 
ral partners such as p and a^. A thorough understanding of 
the dynamics of these mesons in QCD is therefore essential 
for a complete picture of the QCD phase structure. 

In this work we present first results on the properties of 
the chiral partners p and Uj as they emerge from quark- 
gluon fluctuations at high energies. To this end, we study 
dynamical QCD based on [12]. It utilizes the functional 
renormalization group (FRG) approach to QCD [13-22], ex¬ 
tended by dynamical hadronization [15, 23-25]. This way, 
the properties of the hadrons are determined by the underly¬ 
ing dynamics of microscopic QCD and we can conveniently 
describe the transition from quarks and gluons to hadrons 
non-perturbatively without suffering from a fine-tuning of 
model parameters. 

We concentrate on two-flavor QCD at vanishing tempera¬ 
ture and density in Euclidean spacetime and develop a scale 
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dependent effective action that captures the relevant dynam¬ 
ics in both, the quark-gluon phase and the hadron phase, 
on a qualitative level. Furthermore, we extend the dynami¬ 
cal hadronization technique to include vector mesons. The 
present analysis will serve as an starting point for qualitative 
and quantitative in-medium studies of vector mesons. It has 
been demonstrated in [26] for QCD in the quenched limit, 
that quantitative precision is indeed feasible with the FRG 
approach to QGD. 

Even though we work in the vacuum, the renormaliza¬ 
tion group (RG) scale dependence of parameters such as 
the masses reflects their finite temperature behavior. In 
particular, there is a critical scale which separates the 
phases with broken and restored chiral symmetry. This al¬ 
lows us to study the behavior of hadronic parameters as 
they approach the scale of chiral symmetry restoration and 
clarify how mesons decouple from the physical spectrum 
at high energies. The scaling of the low-energy parame¬ 
ters is uniquely fixed from microscopic QGD. The reason 
is that the running of the hadronic parameters is governed 
by an infrared-attractive fixed point as long as the gauge 
coupling is small [24] . This guarantees that the memory of 
the initial conditions of the RG flows of these parameters, 
initiated at large, perturbatively accessible energy scales, is 
lost and the hadronic phase is uniquely determined within 
our truncation. 

By exploiting this fact, we can analyze the validity of vec¬ 
tor meson dominance (VMD) [27]. The idea of VMD is to 
promote the SU{2')i xS!7(2)jj flavor symmetry to a gauge 
symmetry. This way, p and Uj naturally appear as gauge 
bosons [28]. The main advantage is that VMD significantly 
reduces the number of different interactions involving vec¬ 
tor mesons. The obvious conceptual shortcoming is that 
chiral symmetry is only a global symmetry in QCD. Eurther- 
more, even though VMD has lead to accurate predictions 
in some cases at low energies [29], it gives e.g. the wrong 
phenomenology of p and mesons [30]. We therefore 
keep chiral symmetry global and compare our results to the 
corresponding VMD predictions. 

The paper is organized as follows: We motivate our 
ansatz for the quantum effective action used in this work 
in Sec. II. In Sec. Ill we first introduce the ERG and dy¬ 
namical hadronization in the presence of vector mesons. 
We continue with a discussion of the implications of tt— 



2 


mixing within our approach. At the end of this section, an 
outline of our treatment of the gauge sector of QCD is given. 
Our results are presented in Sec. IV. After discussing the 
initial conditions of the RG flows of our model, we present 
our results on the meson and quark masses with particular 
emphasis on the scaling of the p mass towards the chiral 
symmetry breaking scale. Then, we show in more detail 
how the mesons disappear from the physical spectrum at 
large scales. Lastly, we discuss the validity of VMD. We end 
with a conclusion in Sec. V and provide some details about 
the RG flow equations in Sec. A. 

II. THE SCALE DEPENDENT EFEECTIVE ACTION 

We aim at describing the chiral dynamics of two-flavor 
QCD at vanishing temperature and density and in particular 
the vacuum behavior of the vector meson chiral partners p 
and Uj. Furthermore, we want to capture the dynamics of 
QCD over a wide range of scales, taking into account both 
quark-gluon dynamics in the quark-gluon phase at large 
energies and hadron dynamics in the low-energy hadronic 
regime. Thus, we base the construction of the effective ac¬ 
tion we use in this work on well-established renormalization 
group arguments. The starting point of our construction is 
the microscopic gauge fixed action of QCD. Owing to the 
increasing strength of the strong coupling towards lower 
scales, effective four quark interactions are dynamically gen¬ 
erated and become physically relevant. They drive the chiral 
phase transition, as poles in these quark-antiquark scatter¬ 
ing channels signal the formation of bound states and the 
breaking of chiral symmetry. 

The tensor structure of the four-quark interaction chan¬ 
nels is directly related to the quantum numbers of the corre¬ 
sponding mesons that are formed in the respective scattering 
channel. Therefore, we need to include at least those in¬ 
teraction channels, that carry the quantum numbers of the 
mesons we are interested in. In the present case, these are 
the Lorentz-scalar-pseudoscalar iso-scalar-vecor and the 
Lorentz-vector-axialvector iso-vector-vector channels 


and Ayjt- On the mesonic side, these channels correspond to 
the sigma and the pion, and the rho and the respectively. 

The scalar-pseudoscalar channel is the dominant channel 
among all possible quark-antiquark scattering channels in 
vacuum. This has been explicitly checked by considering a 
complete basis of four-quark interactions [26]. This implies 
in particular, that the pions and the sigma mesons domi¬ 
nate the dynamics in the hadronic phase. In this work we 
demonstrate explicitly on the example of vector mesons, 
that there is an emergent scale hierarchy where only the 
lightest mesons, i.e. pions and sigma, can contribute to the 
dynamics of the system at low energies. Thus, the prop¬ 
erties of the heavier meson states in Euclidean space are 
completely fixed by quark-gluon dynamics at large energies 
and pion-sigma dynamics at low energy scales. 

To properly take into account the dynamics in the 
hadronic phase, we model this sector by an effective me¬ 
son potential which in principle includes arbitrary orders 
of mesonic self-interactions. Furthermore, we consider mo¬ 
mentum dependent propagators of the quarks and mesons, 
based on a small-momentum expansion, by including scale 
dependent wave function renormalizations Z^. To account 
for non-vanishing current quark masses, a source term —co¬ 
in the meson sector explicitly breaks chiral symmetry. It is 
directly related to finite current quark masses. As a con¬ 
sequence, pions are massive rather than Goldstone bosons 
and the chiral transition is a crossover. 

To connect the the quark sector with the meson sector, we 
include scalar channel and vector channel Yukawa couplings 
and hyj^. In order to consistently account for the dy¬ 
namical change of degrees of freedom from the quark-gluon 
phase to the hadronic phase, we use dynamical hadroniza- 
tion as it was put forward in [12]. We will elaborate on 
this in the next section. As we will demonstrate there, it 
is inevitable to use this formulation here, since the the 
elimination of the ti —mixing results in manifestly scale 
dependent fields. 

In summary, we use the following scale dependent effec¬ 
tive action: 




( 1 ) 


A A 

+ [(.qqf - (qrsTq)^] - + (qr^rsTg)^] + [qCLs^ff+icr)q] -i- hyj, [qCr^fp'^ -i- u^q] 


+ \zs,k^d^^f + V^)|, 


with the covariant derivative D^ = d^ — iZ ^^where 
is the strong coupling and Z^^^ the gluonic 


wave function renormalization. With this definition the 
covariant derivative is renormalization group invariant. 
We use Hermitian gamma matrices so that [y^, = 25^^!. 
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The commutator for the SU{NJ generators reads [t°, t**] = 
^fabcf-c hence, the trace is positive, Trt“t*’ = f 
are the Pauli matrices. For the field strength tensor we use 
the relation 


We use an 0(4) representation for the meson fields: 

V^ = P‘^f + a^f\ (4) 


F^v= [D^,D^]. (2) 

^A.k Sk 

We work in Landau gauge here, so ^ = 0. For more details 
on the gauge part of our truncation see Sec. Ill C. 

The first line of (1) contains the microscopic gauge fixed 
action of QCD. As mentioned above, we introduced a run¬ 
ning quark wave function renormalization to capture 
some non-trivial momentum dependence of the quark prop¬ 
agator. Aifgiue stands for the fluctuation-induced part of 
the full momentum dependence of ghost and gluon propa¬ 
gators as well as non-trivial ghost-gluon, three-gluon and 
four-gluon vertex corrections, for details see [31-33] and 
Sec. Ill C. 

The four-quark interaction channels and the correspond¬ 
ing Yukawa interactions are in the second line of (1). The 
Yukawa sector arises from the bosonization of the quark 
sector. With dynamical hadronization as explained in the 
next section, these interactions will basically carry the quark 
self-interactions in the quark-gluon regime (see also (A17) 
and (A19)). 

The third line of (1) contains the meson sector of our trun¬ 
cation. With the running wave function renormalizations 
Zgjt and Zy^k for rhn scalar and vector mesons respectively, 
we capture the major part of the momentum dependence 
of the full meson propagators [34] . Furthermore, as we 
explicitly demonstrate in Sec. IV C, the wave function renor¬ 
malizations play a crucial role for the decoupling of the 
mesons at high energies. They are therefore indispensable 
for the identification of the physical meson masses. 

The meson interactions are stored in [/^^((p, V[,), which 
reads 

Uki^,V^) = 

\’T^lk - ‘Fo) + - Vof - 

-^g5,/ctrV^K[V^,Vv]- 

Since within this work our focus is on the qualitative phys¬ 
ical picture of vector mesons in vacuum, we restrict to 
relevant and marginal interactions along the lines of [30] . 
The effect of irrelevant operators, which are potentially non- 
negligible if one is interested in quantitative precision [35], 
is considered elsewhere. Note that I/^ is not an effective 
potential, since it also contains the term ~ with an 
explicit derivative of the scalar meson field. 


and the vacuum expectation value of the mesons is given by 
(po = (0, (JQ^ky ■ We define the so(4) matrices 


(Wjk = 




(T=)=(“|f 



(5) 


with i,j, k e {1,2,3} and ej = (5^;, 52,-, Sj;). They obey the 
following commutation relations: 

[f’o Tj] = 


[T^, T^] — (6) 

[Ti,T/]=ie;,,T,^ 

and therefore = |(T; —T^^) and Tf- = |(T;-I-T;^) form rep¬ 
resentations ofSU[2)i and 51/(2)^}. The mesons transform 
under these flavor rotations as 


v^^uv^uy (7) 

with U = where a and p are arbitrary vectors. 

For more details on the meson sector of our truncation, we 
refer to [30]. 

With this definitions the scalar-vector and vector-vector 
meson interactions in (3) can be written as follows: 

-iguy^f ■ df.'-p = gi,fc X n) ■ 5^71 - cruf • 5^7? 

^^^trV,V,= ^ (Tf^ + a^) ((p'^)^ + (at)^) , 
-^g4,;ttt5^V^[V^,Vj = g 4 _fc [5^p’'-(p''xp’' -Fdf xd]') 

-^g5,)ctrV^V^[V^,Vj = ^ [(p'^xp’'-Fdf xd]')^ 

+ (p'^xdl'-dfxp’')^] . 

( 8 ) 

As it is evident from (3) and mentioned in the introduc¬ 
tion, we do not assume VMD. This gauge principle would 
lead to to the following relations between the different cou¬ 
plings of our truncation: 

glk = g 2 ,k = glk = g5,k and g 3 _;, = 0. (9) 
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By inspection of the renormalization group flow of these 
couplings, we will show that VMD would lead to an over¬ 
simplification of the dynamics of the system. Nonetheless, 
VMD turns out to be a good approximation at low energies, 
see Sec IV D. 

In the present setup, the masses of the quarks and the 
mesons are given by 


K,k = ’T^lk + Kk(^lk’ ( 10 ) 

K,k=^l.k+S3,k(^lk’ 

^i,k=K,k+<^S2,k + 83.k)(^lk- 


We see that the n and the a meson as well as the p and 
Qi meson have degenerate masses in the chirally s)mimetric 
phase which is characterized by aQ ^. = 0. When chiral 
symmetry is broken, this degeneracy is lifted. The mass¬ 
splitting of the scalar mesons is then determined by the 
quartic scalar meson coupling A 4 The mass-splitting of 
the vector mesons is determined by the strength of the 
interaction g 2 jt. Note that, owing to the symmetry breaking 
source c > 0, we are not in the chiral limit. Thus, the chiral 
order parameter CTq./c always nonzero. 

Even though the masses we extract here are the curvature 
masses, it was shown in [34] on the example of the pion 
mass in a quark meson model, that the curvature mass of 
the mesons is almost identical to the pole mass for trunca¬ 
tions that include running wave function renormalizations. 
Thus, as mentioned above, we capture the major part of the 
momentum dependence of the full meson propagators by 
including and and the masses are very close to the 
physical masses. 

We note that even though the action contains massive 
vector bosons, it is not necessary to use the Stueckelberg 
formalism to ensure renormalizability [36]. UV regularity is 
always guaranteed for the functional renormalization group, 
as long as the scale derivative of the regulator decays fast 
enough for momenta much larger than the cutoff scale. 


III. FLUCTUATIONS AND THE TRANSITION FROM QUARKS 
TO MESONS 


A. Functional renormalization group and dynamical 
hadronization in the presence of vector mesons 


Here, we follow the discussion given in [12]. In addition, 
since this work constitutes the first ERG study of vector 
mesons in QCD, we will discuss the implication for the flow 
equations and dynamical hadronization in this case. 

The starting point of the functional renormalization group 
is the scale-dependent effective action at a UV-cutoff 
scale A. In the case of first-principle QCD, A is a large, 
perturbative energy scale and correspondingly F^ is the mi¬ 
croscopic QCD action with the strong coupling constant and 
the current quark masses as the only free parameters. Quan¬ 
tum fluctuations are successively included by integrating out 
momentum shells down to the RG-scale k. This }4elds the 
scale-dependent effective action F;^, which includes fluctua¬ 
tions from momentum modes with momenta larger than k. 
By lowering k we resolve the macroscopic properties of the 
system and eventually arrive at the full quantum effective 
action F = F;t=o- The RG-evolution of the scale-dependent 
effective action is given by the Wetterich equation [37]. 

As we have discussed above, a formulation of the effective 
action in terms of local composite fields is more efficient in 
the hadronic phase of QCD. In order to dynamically con¬ 
nect this regime with the ultraviolet regime of QCD, where 
quarks and gluons are the dynamical fields, we use dy¬ 
namical hadronization as it was put forward in [12]. This 
implies that the meson fields in (1) are RG-scale dependent. 
This yields a modified Wetterich equation, which reads with 
4> = (A, q, q, c, c, ti, ct, p, aj) in a shorthand notation: 


dt^kW = -Tr 


(Ff[4.]+Rf) '-d. 


Rl 


5ct>i 


■ 3t4>i, 
( 11 ) 


where <p = (n,a,p,a^') summarizes the meson fields, is 
the total derivative with respect to the RG-time t = ln(fc/A) 
and the traces sum over discrete and continuous indices of 
the fields, including momenta and species of fields. This also 
includes the characteristic minus sign and a factor of 2 for 
fermions. F®[4>] denotes the second functional derivative 
of the effective action with respect to all combinations of 
the fields. is the regulator function for the field 4>. It 
is diagonal in field space. Note that in order not to break 
chiral symmetry explicitly by our regularization scheme, 
we introduced the same regulators for the scalar mesons 
and the vector mesons respectively. For details we refer to 
App. A. The flow equation can be written schematically as 


In this work we are interested in the dynamical transition 
from UV to IR degrees of freedom. To achieve this, we 
include quantum fluctuations by means of the functional 
renormalization group. For QCD related reviews see [13- 
21]. Furthermore, in order to consistently describe the 
dynamical change of degrees of freedom, we use dynamical 
hadronization [12, 15, 23-25]. This allows for a unified 
description of the interplay between different degrees of 
freedom at different scales in terms of a single effective 
action. 



Here, the first three diagrams represent the gluon, ghost 
and quark contributions to the flow equation and the fourth 
diagram depicts the meson contributions. The black dot 
indicates that the propagators are the full propagators. The 
crossed circles represent the regulator insertion. By inserting 
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the truncation ( 1 ) into the flow equation ( 11 ), one gets a 
closed set of fully coupled RG-flow equations for the scale- 
dependent parameters of the truncation. 

■ di-(pi stands for the modifications of the flow equa¬ 
tion due to dynamical hadronization. The general idea 
is to store four-quark interactions entirely in the meson 
sector at every scale k, i.e. we perform a bosonization at 
every scale [15, 23-25], The hadronization fields cp, which 
carry the quantum numbers of the mesons, become scale 
dependent and can be viewed as hybrid fields; while they 
encode the quark dynamics at large energies, they behave as 
mesons at low energies. This way, a dynamical connection 
between these two distinct sectors is established and the 
transition scale is an emergent scale which is fixed by the 
fluctuations of fundamental QCD. This is in contrast to con¬ 
ventional bosonization by means of a Hubbard-Stratonovich 
transformation, where four-quark interactions are stored in 
the meson sector at a fixed scale and therefore four quark- 
interactions are neglected/mis-counted at other scales [17]. 

The gauge sector as well as the quarks are not affected 
by the hadronization. The scale dependence of the mesons 
is given by their flow df-cp^, which reads for the individual 
mesons 


d,n=A^qY^Tq, 


d,(j=A^qiq, 
5tP''=Bfcqr^iTq, 
d,a^=Bj,qj^Y^zq - Ckd^n. 


(13) 


Note that the structure of the mesons as quark-antiquark 
bilinears becomes apparent in this formulation. A^. and 
Bj. are the hadronization functions. Their precise form is 
not determined a priori. In accordance with the discussion 
above, we fix them such that the fermionic self-interactions 
that drive chiral symmetry breaking and reflect the meson 
content of our theory are stored in the meson sector. Thus, 
the the four-quark interactions are completely absorbed into 
the meson sector, enforcing 


In addition to the quark-bilinear term with the quantum 
numbers of the corresponding meson in (13), the flow of 
the Uj-meson has an additional contribution proportional 
to d^7t. This term arises because the so-called tt—U j mixing 
leads to an additional scale dependence of the meson, 
which has to be taken into account and Axes Q. We will 
elaborate on this point in the next section. 


B. n—Qi mixing 


Spontaneous chiral symmetry breaking leads to a non¬ 
vanishing vacuum expectation value o-q of the cr meson and 
the resulting mixing term 


r 


nai 


gi,k 
J X 


(16) 


implies an off-diagonal two-point function This is 

referred to as n — Uj mixing. Here, we will eliminate this 
mixing by a redefinition of the held, 


-T + 


gl.fcCr 


ml.+[g2,k + g3.k)c^^ 


(17) 


This redefinition of the held renders it explicitly RG-scale 
dependent, 7 ^ 0. Before we discuss the implications of 
this scale dependence, we turn toward the resulting mod¬ 
ifications of the effective action (1). If we plug (17) into 
the truncation ( 1 ), the part of the action leading to the 
mixing term (16) is canceled and various new terms ap¬ 
pear. Since the replacement (17) introduces terms ~ 
the interactions of our original ansatz ( 1 ) receive modifi¬ 
cations with explicit momentum dependence. Within this 
work, we define all running coupling at vanishing external 
momentum, see App. A. Thus, for interactions that are not 
explicitly momentum dependent in our original action, these 
modifications simply drop out of the beta functions. Only 
the meson anomalous dimension and the scalar-scalar 
vector interaction receive non-vanishing modifications. 
The new term ~ {.d^nY yields for the pion wave function 
renormalization 


^t^s,k = 0 and d^?ivj, = 0. (14) 

Note that this formulation eliminates all double- and/or 
mis-counting problems, which potentially occur in models 
including both quark and hadron degrees of freedom [ 12 ]. 
(14) yields the following hadronization functions: 




^k — ~7ri ^tl^^V.k J 


2h 


V,k 


(15) 


where , denotes the scale derivative with fixed hadroniza- 
tion fields. These hadronization functions give rise to modi¬ 
fied running couplings of (1). They are given in App. A. 


^n,k — Zq],- 


glk<k 


S,k ' 


(18) 


While the wave function renormalizations do not enter RG- 
invariant beta functions, their the anomalous dimension 
does. Thus, (18) yields a modified pion 
anomalous dimension. 

The other relevant modification affects the ptitr vertex, 
which now reads 


p(3) 

pnn 


— gl,k 


f 

1 - 
V 


g2,k(^0.k 


m 




(19) 


Since we define the coupling g^. via this vertex, this has to 
be taken into account in the corresponding beta function, 
see App. A. 
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The elimination of the tt —Oj mixing entails a shift of 
the field which includes running couplings. As a con¬ 
sequence, the Qi field becomes RG-scale dependent itself. 
As we have discussed in the previous section, we use the 
dynamical hadronization technique which implies that all 
meson fields are scale dependent. The scale dependence 
of Uj induced by (17) is additional to the one induced by 
dynamical hadronization. The total scale dependence of 
is now given by the RG flow 

= B^qy^Ys'^q - . ( 20 ) 

The first term stems from dynamical hadronization and 
reflects the quark-bilinear nature of the Uj meson. The 
hadronization function is given in Eq. (15). The second 
term is a result of the diagonalization of the meson two- 
point function and according to (17) Q is given by 



Figure 1: The running of the different strong couplings in compari¬ 
son to the 1-loop mnning. 


gi.k^ 


ml + [g2,k + g3,k)o-'^ 


( 21 ) 


In summary, chiral symmetry breaking leads to an off- 
diagonal meson two-point function. Diagonalization leads 
to modifications of the pion anomalous dimension and the 
pnn interaction and introduces an additional scale depen¬ 
dence to the Uj meson. 


regime at scales k <2 GeV, while it leaves the gauge cou¬ 
plings unchanged in the perturbative regime. For more 
details we refer to [12]. 

We also include four-gluon interactions via the four gluon 
vertex g,44 j.. For the sake of simplicity, we choose a semi- 
perturbative ansatz that ensures the correct perturbative 
running and set 


C. Gauge sector 


gA\k — gA^,k ■ 


( 22 ) 


Here, we briefly discuss the gauge sector of our truncation 
(1). We follow the approach to dynamical QGD which was 
put forward in [12] and refer to this work for further details. 

The gauge couplings induced by three-point functions 
play a dominant role in the description of interactions. Most 
importantly, effective four-quark interactions are generated 
by strong quark-gluon interactions during the RG flow in the 
quark-gluon regime. As discussed above, these four-quark 
interactions drive the chiral phase transition and the forma¬ 
tion of mesons in the corresponding quark-antiquark scatter¬ 
ing channels. We therefore solve the full flow equations for 
all three-point functions in QGD, i.e. the quark-gluon ver¬ 
tex gqAq^k’ the three-gluon vertex and the ghost gluon 
vertex 

In order to keep the description as simple as possible, 
despite the rather large number of running couplings we 
consider here, we restrict our analysis to vanishing exter¬ 
nal momentum in the n-point functions. Together with an 
appropriate choice of the regulator function this yields 
analytical flow equations, see App. A. Furthermore, only 
the classical tensor structures are taken into account in 
the gauge sector. It has been shown in [26] for quenched 
QGD, that both non-trivial momentum dependencies and 
non-classical tensor structures in the gauge sector lead to 
important quantitative corrections resulting in larger gauge 
couplings in favor of stronger chiral symmetry breaking. To 
phenomenologically account for these effects, we introduce 
an IR-strength function g(fc) which smoothly increases the 
gauge couplings by a factor ~ 1.17 in the non-perturbative 


This approximation is valid for semi-perturbative scales and 
larger, fc > 1.5 GeV. At smaller scales, non-perturbative 
effects potentially lead to a different running. 

In summary, we consider the following vertices in the 
gauge sector: 


_ y 2 7 c(3) 

^ k ~ ^A,k^q,k gqAq,k ^qAq ' 


^f'’=ZlkgA^k S, 


(3) 
A> ’ 

pG^^) _ y2 2 cC4) 

k ~ ^A,k gA^,k ’ 


pi(cAc) _ y 2 y c(3) 

k ~ ^A,k^c,k gcAc,k ^cAc ‘ 


(23) 


The classical tensor structures ^ are obtained from the 
gauge fixed classical QGD action by 


cM 


54>i... 


(24) 


where we have omitted indices for clarity. The definitions 
of the flow equations of the gauge couplings can be found 
in App. A. The different running gauge couplings follow 
the perturbative running for scales fc > 4 GeV, while non- 
perturbative effects lead to a different running at smaller 
scales. Fig. 1 shows our results for the running gauge cou¬ 
plings. 
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For the gluon and ghost propagators we use the two-point 
functions 


p(2),YM 
^ A/c,fc 


(P) 




(25) 


computed in [38, 39] for pure YM theory as input, and in¬ 
clude matter back-reactions in order to describe unquenched 
QCD [12] . We make use of the fact that within our construc¬ 
tion the propagators enter the flow equation only via the 
corresponding anomalous dimensions 


'nA/c,k 


^t^A/c,k 

^A/c,k 


(26) 


where we only consider external momenta p = k and define 
^A/c,k = ^A/c.ki.^^)- The anomalous dimensions of YM are 
functions of the YM gauge couplings. Thus, in our case, 
where we identified the four-gluon running coupling with 
the three-gluon running coupling, we have 


= ^'^c.k {s™.k’ ) • (27) 

The gluon anomalous dimensions of full QCD consists 
of a pure gauge part rjgiue.it and the vacuum polarization 
induced by quark fluctuations. 


VAk = Pglue.fc + Atja,)c • (28) 

We include the full vacuum polarization at vanishing exter¬ 
nal momentum. For details see App. A. For the pure glue 
contribution of dynamical QCD we make use of (27): Since 
we know that the anomalous dimensions can be expressed 
as functions of the gauge couplings, we replace the the pure 
YM gauge couplings in (27) by those of dynamical QCD, 

Pglue.fc VaA {^ScAc,k^ SA^,k'^ * (29) 

In practice, we use the approximation g™^ = g™^, 
which is valid down to the semi-perturbative regime, 
and parametrize q™ by only one coupling, i.e. q™ 
q™(g™^). Then, we can use a numerical fit function 
q™(g) extracted from the input (27) and insert the cor¬ 
responding QCD coupling. This yields our final expression 
for the gluon anomalous dimension of QCD: 


^Ak = V™igcAc,k) + ^VAk ■ (30) 

We proceed analogously for the ghost anomalous dimen¬ 
sion. There, no direct quark corrections arise in QCD and 
we only need to adapt the input anomalous dimension from 
YM to QCD. As it turns out, r]™ is very well described by 
a function linear in g|. We therefore find for the ghost 
anomalous dimension of dynamical QCD: 


§cAc,k A 

^g™k j 


2 


n™ 

'Ic.k 


(31) 



Figure 2: The RG-flows of the scalar and vector yukawa couplings 
hs k, hyk- The thick and thin lines correspond to different initial 
values of the couplings at different initial scales. 


IV RESULTS 

A. Initial Conditions 


We initiate the RG flow of the effective action (1) at 
the initial scale A = 20 GeV and therefore deep in the per- 
turbatively accessible quark-gluon regime. Hadronic de¬ 
grees of freedom will become relevant at much lower, non- 
perturbative scales. Owing to the dynamical hadronization 
procedure, the RG flows of our model are completely fixed 
by specifying the free parameters of the microscopic gauge 
fixed action of QCD, i.e. the strong coupling and the cur¬ 
rent quark mass. Even though we choose a priori different 
flow equations for the strong couplings a^s, the 
gauge principle enforces them to be identical in the per¬ 
turbative regime. The initial value of the strong couplings 
implicitly sets the scale, and we choose a^ l^ = 0.163 for 
i = qAq,A^, cAc, which corresponds to A 20 GeV. Indeed, 
Fig. 1 shows, that the gauge couplings are all identical and 
follow the perturbative 1-loop running down to fcf%i4GeV, 
where non-perturbative effects start to induce differences 
between these couplings. 

The current quark mass is related to the explicit 
symmetry breaking parameter c in (3) via 


m 


uv 

q 



(32) 


We choose for the renormalized parameter = 3.9 GeV^, 

which yields an pion mass in the IR of M„ q = 137.5 MeV. 

We note that the physical parameters are rescaled with 
appropriate powers on the wave function renormalizations 
to ensure RG invariance, see App. A and in particular (A3) . 
The physical (or renormalized) quark and meson masses 
are defined as 


^q,k — 


^q,k 


and 


^cl>,k — 


T-4,,k 


7 I /2 ■ 
"4>,k 


Vc.k = 


(33) 
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k [GeV] k [GeV] 

Figure 3: Renormalized masses (33) as a function of the RG-scale. Masses in the shaded area are larger than the cutoff scale and therefore 
decoupled from the dynamics. 



Figure 4: Effective propagators (A9) of the matter fields as a 
function of the RG scale. They are a measure for the effective 
strength of the fluctuations of the fields. 


With slight abuse of terminology, we refer to rrif^ as bare 
mass. They are given in (10). 

The initial conditions of the mesonic parameters can be 
chosen arbitrarily. In the regime of weak gauge coupling, 
the flows of these couplings are governed by an infrared- 
attractive fixed point [24]. Thus, as long as the initial scale 
is large enough, we find unique solutions for the meson 
parameters at low energies. This is given as long as the 
initial meson masses are chosen larger than the UV-cutoff 
scale, > A. That way, the mesons do not contribute to 
the dynamics of the system at high energies. Furthermore, 
to ensure that our initial conditions correspond to QCD, the 
ratio J^\ip K much smaller than A“^. It 

corresponds to the four-quark coupling ^siv,h ihe initial 
scale. A large initial value of the four-quark coupling would 
describe a gauged Nambu-Jona-Lasinio model with strong 
coupling, rather than QCD. 

The independence of the IR-physics on the initial values 


of the meson sector is demonstrated in Fig. 2 and Fig. 7. 
There, we have chosen initial values at different initial scales 
(10 and 20 GeV) that differ by many orders of magnitude 
and one nicely sees that the initially different trajectories are 
attracted towards a unique solution in the hadronic regime. 


B. Masses 

In Fig. 3 we show our results for the quark and meson 
masses. The left figure shows the masses over the full range 
of scales we consider here, while the right figure shows the 
region for k < 600 MeV. For scales k > 400 MeV all mesons 
are decoupled from the flow. At these scales the dynamics 
are driven completely by current quarks and gluons. At 
about 400 MeV, the degeneracy of the n and o masses as 
well as the p and masses is lifted due to chiral symmetry 
breaking, n, o and the constituent quarks are the dynamical 
degrees of freedom in this region. The vector mesons are 
always decoupled. Thus, the vacuum structure of the vector 
mesons is determined by quark and gluon fluctuations at 
large scales and the fluctuations of the lightest mesons, the 
71 and <j, at lower scales. 

This is also shown in Fig. 4, where the effective propa¬ 
gators also reflect this scale hierarchy. They are defined 
in App. A and in particular (A9) and are a measure for 
the strength of the fluctuations of the fields. Vanishing of 
the effective propagator of a field implies that this field 
does not contribute to the dynamics of the system. Thus, 
we see that at large energy scales the quarks are the only 
dynamical matter fields. There is only a relatively small win¬ 
dow, 100 MeV <k< 500 MeV, where the scalar mesons are 
d)mamical and one nicely sees that the sigma mesons decou¬ 
ples earlier that the lighter pions. Vector meson fluctuations 
are always negligible. 

Indeed, an explicit calculation of a complete set of four- 
quark interactions in Euclidean spacetime shows that the 
scalar-pseudoscalar channel is the dominant channel [26]. 
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Figure 5: Scalar and vector wave function renormalization as 
function of the RG scale. We normalized them to be 1 at the IR 
scale fc = 30MeV. 


This implies that the only relevant meson degrees of freedom 
in vacuum are n and u. We note that this picture will 
change in Minkowski space, since different channels will 
become relevant as soon as the momentum is close to the 
corresponding mass pole. 

The behavior of the masses as function of the RG scale k 
reflects their scaling with temperature T. In particular, the 
running of the masses at vanishing temperature is qualita¬ 
tively very similar to the temperature dependence at fc = 0 . 
In Fig. 3 we see that the p mass is almost constant in the 
hadronic phase and only grows slowly close to the pseudo- 
critical scale. This in line with the in-medium behavior of the 
p mass observed within effective field theory studies, which 
reproduce the data on vector meson spectral functions and 
dilepton spectra very well [ 6 ]. For a sensible comparison, 
however, we also need to compute the in-medium modifica¬ 
tions of the masses within our QCD-based approach. 

Our predicted masses for the vector mesons show a quite 
large discrepancy from the observed masses. We find for the 
renormalized p mass Mp q = 990 MeV, which is about 29% 
larger than the observed mass of 770 MeV [40]. For the 
mass we find = 1077MeV, which is about 15% smaller 
than the observed mass of 1260 MeV The value of the p 
mass is fixed mainly by the fluctuations in the quark-gluon 
sector. This can be seen from the definition of the masses, 
( 10 ), the observation that fhyj^ runs only very little in the 
hadronic regime and that very small (see Fig. 7). 

Thus, the strength of the four-quark interaction which 
is determined by the strong coupling, essentially fixes the 
p mass. Furthermore, according to (10), the mass-splitting 
of p and Uj and therefore the mass of Uj is determined by 
the flow of the hadronic sector at low energies. We note 
that the situation is different for n and cr: the mass of the 
pion is fixed by its nature as a (pseudo) Goldstone boson 
and the strength of explicit S 5 mimetry breaking in terms of a 
finite current quark mass. In any case, the mass-splitting of 
chiral partners in the phase with broken chiral symmetry is 
sensitive to the quality of our truncation in the hadronic sec¬ 


1.2 



0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 

k [GeV] 


Figure 6: Bare masses of the mesons, j., see (10) and 

(33). Note that we normalized the wave function renormalizations 
to be 1 in the IR at fc = 20MeV here. 


tor. Thus, the small mass of Uj is a signal for a shortcoming 
of our truncation there and may be related to momentum 
dependencies that were taken into account insufficiently. 
The large value for the p mass can also be attributed to 
the insufficient inclusion of momentum dependence, but in 
the quark-gluon sector and in particular in the four-fermi 
interaction Ryj.- evaluate this interaction at vanishing 
external momentum, see App. A, and it is possible that we 
underestimate its strength this way. A larger leads to 
smaller vector meson masses. Since the present work is the 
first study in this direction, aimed at capturing the qualita¬ 
tive features, we defer a thorough quantitative analysis to 
future work. 


C. Decoupling of the Mesons 

Mesons are not present in the quark-gluon plasma. In a 
formulation of the dynamics of QCD on a very wide range 
of scales in terms of one scale dependent effective action, 
as in the present work, however, hadronic parameters are 
necessarily a part of the action also at very large scales. As 
we have demonstrated here, the meson masses are much 
larger than the cutoff scale in the quark-gluon regime and 
therefore the mesons are completely decoupled in this phase. 
We want to emphasize that this physically desirable picture 
is achieved with dynamical hadronization. The decoupling 
of the mesons is triggered by a rapid fall-off of the meson 
wave function renormalizations Zg/yj. at the pseudocritical 
scale. Their running is shown in Fig. 5. While they stay 
almost constant in the hadronic regime, they rapidly fall-off 
at about 400 MeV The scalar meson wave function renor¬ 
malization Z 5 fc drops about eight orders of magnitude and 
that of the vector mesons, Zyfc, about seven orders of mag¬ 
nitude towards the UV The fastest drop-off is in the vicinity 
of the pseudocritical scale 400 MeV The reason is 
that quark fluctuations decrease the meson wave function 
renormalizations in the quark-gluon regime. Their flows 
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are proportional to the corresponding squared Yukawa cou¬ 
plings, oc -hl^yj^Zsiv,k, at scales k > k^, resulting 

in large negative beta functions, see Fig. 2. Since the wave 
function renormalizations are the coefficients of the kinetic 
terms in the effective action (1), their vanishing implies that 
the mesons become auxiliary fields and are therefore not 
part of the physical spectrum at large energy scales. 

This is reflected in the behavior of the bare masses, i.e. 
the masses without rescaling with the wave function renor¬ 
malizations, (33), shown in Fig. 6. The 

bare masses would not decouple in the quark-gluon regime: 
while they do not differ from the renormalized masses 
(Fig. 3) in the hadronic regime where the wave function 
renormalizations are almost constant and of order one, they 
are constant in the quark-gluon regime. Thus, at large scales 
the bare meson masses are always much smaller than the 
the cutoff scale. Without the rapid fall-off of the meson 
wave function renormalizations, the mesons show no decou¬ 
pling, resulting in an unphysical high-energy phase. Note 
that the constant bare masses imply in particular that the 
running of the physical masses is exclusively driven by the 
anomalous dimensions of the corresponding mesons at large 
energy scales. 

Since the wave function renormalizations only enter the 
set of flow equations through the corresponding anomalous 
dimensions, the flow equations for the wave function renor¬ 
malizations do not need to be integrated for the solution 
of the system and all results are independent of the initial 
values Zj^. For illustration purposes, we have chosen the 
initial conditions such that Z^jy^ = 1. 

D. Vector Meson Dominance 

The principle of vector meson dominance (VMD) entails 
that the SuINj:\ x SU[Nf)y flavor symmetry is treated 
as a gauge symmetry. In this case, the vector and axial- 
vector mesons appear as gauge bosons of the scalar and 
pseudoscalar mesons. This simplifies the effective action in 
the hadronic sector, since the gauge principle significantly 
restricts the number of possible different interactions and 
there is only one running coupling for interactions involving 
vector mesons. Here, we do not apply VMD. As a conse¬ 
quence, we have a priori different running couplings gi-s^k> 
while VMD implies 

Si.k = g4,k = = \f^k and gs = 0. (34) 

As we have discussed above, the advantage of our approach 
is that the hadronic parameters are uniquely determined 
by the dynamics in the quark-gluon phase, i.e. microscopic 
QCD. Thus, even though we have a large parameter space in 
the meson sector, model parameter tuning is not necessary. 
This allows us to study the validity of VMD in an unbiased 
way by comparing our results to (34). In Fig. 7 we show 
our results for the running of gi_ 5 _fc. 

Or results show that, while VMD does not hold exactly, 
it is a good approximation. In particular the couplings 



Figure 7: The RG-flow of the vector-vector and vector-scalar meson 
couplings gi _5 It. The thick and thin lines correspond to different 
initial values of the couplings at different initial scales. Vector 
meson dominance would imply (34). 


g 2 ,k, g 4 ,k and gs are very close together. Only is con¬ 
siderably larger than the other couplings. If we define the 
error one would make by assuming VMD by the standard 
deviation of these couplings in the IR, we find it to be 
about 16% of the mean average of these couplings, g^j^, 
which is explicitly forbidden for local chiral symmetry, is 
well approximated by VMD. It is very close to g^ ^. = 0 at 
the pseudocritical scale and assumes only a small finite 
value at lower scales. The flow of g^^k is proportional to 
the chiral order parameter Uqj^. Thus, with large positive 
anomalous dimensions, the renormalized coupling g^ ^. is 
driven to values very close to zero at large energy scales. 

The construction of our effective action is based on a small 
momentum expansion (derivative expansion) and we define 
all running coupling at vanishing external momentum, see 
App. A. The momentum scale of our results is therefore given 
by k. Thus, our findings in the hadronic regime correspond 
to small momentum scales k < 400 MeV. A comparison 
of effective field theory predictions assuming VMD with 
experimental results for the electromagnetic form factor of 
the pion show that they agree within 10-20% accuracy at 
momentum transfer < IGeV^ [29]. Thus, our results 
for the validity of VMD are in very good agreement with 
phenomenological findings. 

We note again that the thick and thin lines in Fig. 7 
correspond to very different initial conditions for the flow 
of the couplings. The flows in the hadronic phase as well as 
the final value of the couplings in the IR are prediction of 
our analysis without any model parameter fixing. 


V CONCLUSIONS AND OUTLOOK 

A thorough understanding of the dynamics of vector 
mesons in QCD is essential for our understanding of the 
phase structure of strongly interacting matter. Since these 













11 


low-energy degrees of freedom ultimately derive from mi¬ 
croscopic QCD, the dynamical connection between the high- 
and low-energy sector of QCD needs to be captured. To this 
end, we have presented the first functional renormalization 
group study of vector mesons in QCD. Our focus was on 
how the dynamics of the lightest meson chiral partners, rr, 
o and p, Ui, emerge from the dynamics of quarks and glu¬ 
ons. We have developed a scale dependent effective action 
that captures the dynamical transition from the quark-gluon 
regime to the hadronic regime, including vector mesons, in 
a qualitative manner. The key ingredient is the dynamical 
hadronization technique, which allows for a consistent de¬ 
scription of the transition from high-energy to low-energy 
degrees of freedom. This entails in particular that the prop¬ 
erties of the hadronic regime are fixed by the quark-gluon 
fluctuations at high energies. Thus, no fine-tuning of model 
parameters is necessary and e.g. the meson masses and 
the running of the mesonic parameters can be viewed as 
predictions from first-principle QCD. 

We have demonstrated explicitly that, within this Eu¬ 
clidean formulation, there is an intriguing scale hierarchy 
emerging, where the hadronic contributions to the dynamics 
of the heavier mesons are determined solely by the fluctua¬ 
tions of 71 and cr. 

The masses of p and Uj are almost constant and only 
slightly grow towards the pseudocritical scale. Since the 
behavior of the masses as a function of the RG-scale re¬ 
flects their finite temperature behavior, this gives a hint for 
the in-medium scaling of these masses. Our predictions 
are in agreement with the findings of phenomenologically 
motivated effective models. 

We have emphasized the important role that the meson 
wave function renormalizations play for the decoupling of 
the meson degrees of freedom at high energies. They fall-off 
many orders of magnitude in the vicinity of the pseudocriti¬ 
cal scale. This triggers a rapid growth of the renormalized 
meson masses and the mesons become auxiliary fields in 
the quark-gluon phase. 

Since the properties of the mesonic parameters in our 
model are fixed by the QCD flow, we have been able to 
make an unbiased analysis of the validity of vector-meson 
dominance. Our results show that while VMD does not hold 
exactly, it is a good approximation within an accuracy of 
about 16% at small momentum scales. This is in agreement 
with phenomenological findings. 

In this study we focused on qualitative features and given 
the lack of quantitative precision, in particular for the mass¬ 
splitting of the mesons in the chirally broken phase, there 
is a lot of room for improvement and refinement. In partic¬ 
ular the extension of our truncation in the hadronic sector 
and a thorough analysis of momentum dependencies are 
important next steps. 

This work serves as a starting point for the study of the 
in-medium modifications of the vector mesons and their 
spectral functions within functional renormalization group 
methods for QCD. 
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Appendix A: Running couplings 

In this section we provide some details about the RG flow 
equations of the running couplings of our truncation (1) 
together with the modifications that result from dynamical 
hadronization (13). Due to excessive length of the explicit 
flow equation, we only present their definitions here. For 
the derivation of most of the equations we used an extension 
of DoFun [41] which utilizes Form [42] and FormLink [43]. 
This extension was developed and first used by the authors 
of [26]. With the truncation (1), the Wetterich equation 
(11) and the definitions given below, the flow equations of 
the couplings are uniquely specified. 

Physical parameters are RG-invariant quantities. To 
achieve this, all fields are rescaled with their respective 
wave function renormalizations, $ ^ ^7^$, and all cou¬ 
plings are rescaled with appropriate powers of the wave 
function renormalizations accordingly, see below. This en¬ 
tails in particular, that the wave function renormalizations 
enter the flow equations only through the corresponding 
anomalous dimensions. 


The physical, i.e. RG-invariant, parameters of the action are 
defined as 

’ (A2) 

I 7'^c 7^*#’ ] 

\^A,kq,k^c,k <t>M V,k J 

where 

^ {SqAq,k’ SA^,k^ SA'^,k^ ScAc,k> ^S,k> ^V,k^ ^S,k^ 

(A3) 

hv,k> kngj^, myj^, Vj^, gi_5,;t} , 

is one of the running couplings of our truncation (1). 
n^^, rig, n^, n^, ny are the numbers of gluon, quark, ghost, 
scalar meson and vector meson fields respectively, that are 
attached to the coupling 7)c- The physical masses are given 
by (33). Note that the definition of the gluonic vertices (23) 
implies that the gauge couplings are already RG-invariant. 
Thus, in that case = gj. and we omit the bars. 
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We use 4d regulator functions of the form 

Rlip^) = r^Piirpip^/k^), 
KiP^) = Zc,kP^rBip^/k^), 


(11). Thus, only the first term of the propagators in (A8) 
contributes in the final flow equations. We therefore de¬ 
fine the effective propagators relevant for the flows of the 
physical quantities as 


(A4) 




1 

1 + [M^_k/kf ■ 


(A9) 


R^(p^) = Z^j^p^rgip^/k^), 

= ■^v,icP^''B(p^/fc^)n-^, 

with the transversal projection operator 


^ixv 


P|LlP-V 


(A5) 


For the bosonic and fermionic regulator shape functions 
and rp we use the optimized shape functions [44] 


rpM 


i-ijeci-x), 


rpix) 



0 (1-x). 


(A6) 


The flow equations presented in the following are derived 
using these specific regulators. They have the advantage, 
that the loop-momentum integration can be performed ana¬ 
lytically and, consequently, all beta functions can be given 
in analytical form. Furthermore, we work in Landau gauge, 
fix the Euclidean spacetime dimension to d = 4 and color 
and flavor are fixed to iV^ = 3 and Nf=2. 

First, we explain the effective propagators used in 
Sec. rv B and in particular Fig. 4. The propagators in mo¬ 
mentum space are of the form 


Vanishing implies that the field 4> does not contribute 
to the dynamics of the system. 

We proceed with the definition of the flows of the gauge 
couplings. The explicit form of the flow equations is given in 
[12] . Here, we only present our definitions for completeness. 
As we have discussed in Sec. IIIC, we compute all three- 
point functions of QCD, but restrict them to have only the 
classical tensor structure. We therefore define the flow of 
the quark-gluon vertex as 


='!'o 

(AlO) 


^tSqAq,k 


limTr 


8Nj:[N^ - 1 ) P-o 




where 4>o = (0,0,0,0,0,0, aQj^, 0,0) is the vacuum expec¬ 
tation value of the mean field $ = (A,q,q,c,c, n,a,p,a^'). 
The trace runs over all external indices and includes a loop- 
momentum integration. The limit denotes that all external 
momenta are set to zero. We define the three-gluon vertex 
via the projection 


^tSA^.k — 


lim 


12iV,(iV2 - 1) p^o 


(All) 


r 


Tr 


V 


^ nvPcrf 


abc 


5^8,r. 


5A(p)“5A(-p)j5A^^(0)^ 


GB.kCq) = 


Gpr(q) = 


ZB.kq^ (1 + rgiq^/k^)) + ml ,^ ’ 


^Ik^^ (1 + '■pCqVfc^)) + mlk 


(A7) 


for bosons and fermions respectively. For the specific choice 
of regulator shape functions (A6), they read 


Gfi./tCq) 


0 (fc 2 -q 2 ) 0 (q 2 -fc 2 ) 

Zg^kk^ + mlj^ Zg^^q^ + mlj^ ’ 


Gp.fcCq) — 


0 (fc^-q^) 0 (q^-fc^) 


(A8) 


To be consistent with the low-momentum expansion our 
construction of the effective action is based on, we define all 
running couplings at vanishing external momentum. This 
entails that all integrands of the loop-momentum integra¬ 
tions are proportional to 0(fc^ — q^) which stems from the 
scale derivative of the regulator in the flow equation 


Since as an approximation we evaluate all flow equations 
at vanishing external momentum, the ghost-gluon vertex 
ScAc,k only has canonical running. The diagrams that con¬ 
tribute to the beta function are proportional to the external 
momentum and therefore vanish here and we are left with 

1 . . 

StScAc = - [VA,k + 2Pc,k) gcAc- (A12) 

Since we approximate the four-gluon vertex with the three- 
gluon vertex, see (22), we do not need a separate equation 
for this coupling. 

Next, we discuss the flow of the four-quark cou¬ 
plings. Here, we consider two channels, the scalar- 
pseudoscalar channel with coupling and the iso-vector- 
iso-axialvector channel with coupling Some caution 
is advised when four-fermion interactions are included in 
the effective action. A specific quark-antiquark interaction 
channel can always be expressed as an linear combination 
of different interaction channels with two spinor fields in¬ 
terchanged. This can potentially lead to ambiguities in the 
corresponding bosonized models since different sets of com¬ 
posite states can be related to one and the same fermionic 
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action (see e.g. [20]). This is known as the Fierz ambigu¬ 
ity. While this ambiguity can lead to large uncertainties in 
mean-field calculations, appropriate approximations that 
go beyond mean-field in RG studies can minimize these 
uncertainties [45]. 

Indeed, as explicit calculations considering the RG flows 
of a Fierz-complete basis of four-quark interactions have 
shown [26], the scalar-pseudoscalar (S—P) channel is the 
dominant channel in vacuum, while all other channels are 
strongly suppressed compared to this channel. Furthermore, 
the dynamical hadronization of only the (S —P) is sufficient 
to render all four-quark interaction channels finite at the 
chiral phase transition. Thus, the error we make from not 
using a Fierz-complete basis is expected to be small. We 
therefore restrict our model to contain only two physically 
relevant channels. In order to study the properties of the 
corresponding composite fields, we dynamically hadronize 
both channels here. 

We define the running coupling of the scalar-pseudoscalar 
channel via the projection 

I 1 

^tL^s)c=-;-1 (A13) 

8NfN^i2NfN, + l) 


X lim Tr 

p^O 




543 ,r, 




4.=4.„ 


where A, B, C, D abbreviate the color, flavor and spinor 
indices of the quarks. For the vector-axialvector channel we 
choose 


, ^V,k — ~ o 


I ^S,k 


(A14) 


- lim Tr 

p^O 



545 ,r, 






Taking dynamical hadronization into account, the total flow 
of the renormalized scalar Yukawa coupling is 

dthk = d,\^hs,k - , (A17) 

where ^ndA^ is given by (15). According 

to (10), the scalar channel Yukawa coupling defines the 
quark mass. 

We define the vector Yukawa coupling hyj^ via the pqq 
three-point function as 


(A18) 

where contractions over the remaining indices with Kro- 
necker deltas is implied. Splitting the flow into the con¬ 
tributions with and without dynamical hadronization, we 
find 

3thk hk - k-^Mlk^k, (A19) 

with and B^. given by (15). 

We want to emphasize that the modifications of the 
Yukawa couplings in (A17) and (A17) proportional to Aj. 
and B^. are crucial for the dynamical hadronization proce¬ 
dure. They guarantee that the ratio h^/v^k^^%p k replaces 
the four-quark interactions which vanish due to dy¬ 

namical hadronization, in the quark-gluon phase. This way, 
the modified Yukawa couplings capture the relevant quark- 
gluon dynamics at large energy scales, while they act as the 
usual Yukawa couplings in the hadronic regime. 

Next, we discuss the mesonic couplings of our truncation. 
They are not modified by dynamical hadronization. We 
define the running of the chiral order parameter aQj^ via 
the pion two-point function as 


dthv,k = 


16N,[N^ ■ 


— limTr r,,T— --r I 

1 ) p->o ^ 5pf^5q5q J 


5^9,r, \ 


with the projection operator 


jpjABCD 


1 

l92NfN, 


^ab^cdT 


O-aO-b 


yCtcCtp, 

I /i 


(A15) 


Here, c,n is the spinor index of the respective quark 
field. The Kronecker deltas are summed over the remaining 
color and flavor indices. 

We note that these projections give the flow equations 
for scale-independent meson fields, i.e. without dynami¬ 
cal hadronization. Dynamical hadronization enforces (14). 
Nevertheless, the flows of the four-quark interactions de¬ 
fined in (A13) and (A14) play a major role for the dynamics 
of the system and enter the hadronized flow equations in 
the meson sector via the hadronization functions (15). 

Following the discussion in [35], we define the scalar 
Yukawa coupling hgj^ via the quark-antiquark two-point 
function as: 


-i f 

9ths k =-1™ Tr 5jyn -— 

ANfN, cjq p-o ^ Sq^Sq^ 


(A16) 


9t<^o,k = - [ VfcCr Q - — 


lim Tr 5, 


■ 

f 


1 p->o 


5^9,r, 

SKiSTlj 


(A20) 


with the adjoint flavor indices i,j. The flow of the scalar 
four-point function is defined as follows: 


9tVk 


- 1 p-o 


limTr 5,,5 


'ij^kl 


5^9,r, A 


(A21) 


with the adjoint flavor indices i,j, k,l. 

The explicit symmetry breaking term c is a source term 
and therefore drops out of the flow equation. The RG- 
invariant coupling therefore only runs canonically, 

1 

9t^k = 2^S,k^k ■ 




(A22) 
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The meson masses are defined as the momentum inde¬ 
pendent part of the corresponding two-point functions. For 
the scalar mesons, we need the flow of which is given 
by 


and the projection operator for 

pf' = —j+ 5ik5u] ■ (A28) 


dtml,= 


S,k j^2 , 


— limTr | 5, 

1 p-*o 


5^ 

' SrciSTij 


(A23) 


We cannot define the flow of the vector meson mass pa¬ 
rameter TUyj^ independently of other couplings, since we 
have to either project on the p or the mass, which gives 
contributions from other couplings in the chirally broken 
phase according to (10). We choose to project on the p 
mass and find 


The vector meson self-interactions and g^ j^ are de¬ 

fined as 


3tg4,k = 


— I 


lim ■ 


6NfiN^ - 1) P^o dp 

f 

^ Pa^jdy^ijk 


(A29) 


5^8,r. 


V 


5p“(p)5p; (-p)5p[(0) 






limTr 


- 1 ) p -6 


S^d.Vk 

^ ^ _ I- ^ 

^ jLi 




M 


(A24) 




The flow Sfgs is defined below in (A26). 

For the definition of the three-point function gyj^ we 
choose the ptiTi vertex and find; 


and 


dtgs.k = 


24A//(A// 


— lim 

1 ) P^O 


(A30) 


X Tr 


f 

^ajS^yS^ij^kl 

V 


5^3,r, 

5p^5p^5pl5pf ^ 


'!’=4’o 


—i . 8^ 

^tgl,k /' at 2 PI 2 

2Nf{Nj - 1) p^o 8p‘^ 

5^8,Y 


Tr 


k 


5pr(-p)57t/p)57rfc(0) 


+ ^l.k^t 


Sl,kS2,k 


- g2,kO-o,kCk ■ (A25) 


V 


As we have discussed in Sec. Ill B, the elimination of the 
71—mixing leads to two types of modifications of the pnn 
vertex. The first stems from the modifications of the action 
due to the replacement (17) and leads to a modification of 
this vertex given by (19). The first term in the third line of 
(A25) cancels the additional term in the flow to ensure that 
we compute the flow of gy ^ and not of (19). The second 
modification stems from the scale dependence of Uj that is 
introduced by (17). This leads to the second term in the 
third line of (A25) which follows from (13). 

We define the couplings g 2 ^k g^^k via the flow 


^tg2/3,k 


limTr 

p^O 


f S'^d.tk 

5 P'-''^*- — - 

^ 5jii57ij5p^5p 



4>='to 

(A26) 


with the projection operator for g 2 ^k 


■pijkl 

Si 


4Nf[Nf + l) 


Nf - 1 


5.. 5, 


^ik^jl 


(A27) 


Finally, we discuss the wave function renormalizations. As 
mentioned before, in a RG-invariant formulation they enter 
the flow equations only via the corresponding anomalous 
dimensions (Al). The ghost anomalous dimension and the 
gauge part of the gluon anomalous dimension are discussed 
in Sec. Ill C. the quark contribution to the gluon anomalous 
dimension i.e. the vacuum polarization, is computed 

from 


^PA,k 



(A31) 


Again, the lines correspond to the full propagators and the 
dots to the full vertices. The explicit equation can be found 
in [12]. 

The quark anomalous dimension is computed from 


-1 


(A32) 


8^ f 5^8,tk A 

^P™ 8 p^'^^ {j^^^ 5 q(j)) 5 q[-p) J 

where contraction of external color, flavor, and spinor in¬ 
dices is understood. 

The scalar meson anomalous dimension rj^ k has to be de¬ 
fined via the pion-pion two-point function. Using the sigma 
meson two-point function leads to additional contributions 
to the flow with couplings that correspond to a higher order 
derivative expansion. We define the anomalous dimension 
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for scale-independent fields as 


Is.k 


<i> 


-1 

2(iv; - l)Zs,t 


(A33) 


X lim ^:-^Tr | 5 
P^o dp^ 


5ni[p)5n:[-p) 


4'='!'n 


It receives modifications from the scale dependence of 
from the elimination of the tr—Q j mixing. The full anomalous 
dimension then is 


Vs.k= Vs.k\^- g 2 ,k o-Q.fcQ (A34) 

where the second term follows from (13). We use the rho 
meson to define the vector meson anomalous dimension 
and find 


We emphasize that due to chiral symmetry the definition 
of the mesonic couplings in terms of n-point functions is not 
unique. We have explicitly checked that different projection 
procedures give the same results as long as they are equiva¬ 
lent by chiral symmetry. However, some caution is advised 
since seemingly equivalent definitions may give different 
results. The reason in those cases is that that inappropriate 
projections may contaminate the flows with contributions 
that are not part of the truncation. For example, a defini¬ 
tion of vj. via the sigma meson four-point function instead 
of (A21) gives additional contributions from diagrams that 
are related to the flow of the 6-meson interaction. Another 
example is which is mentioned above (A33). To find 
appropriate projection procedures one therefore has to keep 
extended truncations, such as general field-dependent cou¬ 
plings, in mind. 


Vv.k = 


-1 

6iNf - l)Zv,fc 


(A35) 


X lim 7 -^Tr 
P^o dp^ 


( 5^d,r. 

- U - — - 

V 5pf(p)5pJ(-p)^ 
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